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Abstract 

The coupHng of gravity to matter is explored in the hnearized gravity hmit. The usual derivation 
of gravity-matter couplings within the quantum-field-theoretic framework is reviewed. A number 
of inconsistencies between this derivation of the couplings, and the known results of tidal effects 
on test particles according to classical general relativity are pointed out. As a step towards re- 
solving these inconsistencies, a General Laboratory Frame fixed on the worldline of an observer 
J> ■ is constructed. In this frame, the dynamics of nonrelativistic test particles in the linearized grav- 

. ity limit is studied, and their Hamiltonian dynamics is derived. It is shown that for stationary 

I metrics this Hamiltonian reduces to the usual Hamiltonian for nonrelativistic particles undergo- 

' ing geodesic motion. For nonstationary metrics with long-wavelength gravitational waves (GWs) 

O . present, it reduces to the Hamiltonian for a nonrelativistic particle undergoing geodesic deviation 

. motion. Arbitrary-wavelength GWs couple to the test particle through a vector-potential-like field 

bX). 

Na, the net result of the tidal forces that the GW induces in the system, namely, a local velocity 



^ ■ field on the system induced by tidal effects as seen by an observer in the general laboratory frame. 

Effective electric and magnetic fields, which are related to the electric and magnetic parts of the 
Weyl tensor, are constructed from Na that obey equations of the same form as Maxwell's equations 
. A gedankin gravitational Aharonov-Bohm-type experiment using Na to measure the interference 
of quantum test particles is presented. 
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I. INTRODUCTION 



At first glance, it would seem that the Hamiltonian dynamics of nonrelativistic, classical 
test particles in the linearized gravity limit has been thoroughly studied, and is well under- 
stood. Indeed, in this limit gravitational waves (GWs) are often treated as simply a spin-2 
gauge field propagating in flat Minkowski spacetime [ll , and the coupling of GWs to matter 
would seem to follow naturally. This determination would be premature; we show in this 
paper that such an approach obfuscates the underlying physics of the system, and overlooks 
the surprising links between gravitational waves, vector potentials, and gauge symmetries. 

Much of our current understanding of the coupling of matter to gravity comes from at- 
tempts at constructing quantum gravity (QG) 

mm 

, and from the theory of quantum 
fields in curved spacetime (QFCS) js, |^. Once the Lagrangians for various elementary 
particles — both gauge and nongauge — were determined in fiat Minkowski spacetime, their 
extension to curved spacetimes was a natural next step. To make this extension of the 
flat spacetime Lagrangians to curved spacetimes, a number of seemingly natural assump- 
tions were typically made then, and are still being made now. The expectation is that the 
experience and intuition gained from constructing quantum field theories (QFTs) in flat 
spacetime will serve as useful guides in constructing QFTs in curved spacetimes. Thus, flat 
spacetime Lagrangians for bosonic fields are promoted to curved spacetimes by replacing 
the Minkowski metric with the metric for a curved spacetime, the partial derivative with 
the covariant derivative, and the Lorentz-invariant integration measure with the general- 
coordinate-transformation-invariant integration measure. The extension of fermionic fields, 
such as spin- 1/2 and -3/2 fermions, follow in much the same way once a tetrad frame is 
chosen. The Hilbert action from classical general relativity (GR) is used for the gravity 
component of the theory, and the metric for the spacetime g^^, is identified with the grav- 



itational fie^ 
chosen — in 



d's degree of freedom. A classical backgroiind_metric g^^ for the spacetime is 

it was either an asymptotically 



l| it was flat Minkowski spacetime, and in 0, 
flat spacetime or a spacetime with finite spatial extent — and the propagating component of 
the gravitational field is extracted from the theory by considering fluctuations h^^ about g'^j^ 
combined with a suitable gauge (coordinate) choice. These fluctuations — representing grav- 
itational waves (GWs) in GR and gravitons in QG — are then expanded about (7°^ assuming 
that h^i, = g^i, — g^ are small compared to (7° , and are subsequently treated as simply 
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another spin-2 non-Abelian gauge field propagating in the background spacetime. By also 
expanding the metric terms about (7°^ in the Lagrangians for the matter fields, one obtains 
terms that couple GWs — or gravitons in QG — with matter. These interaction terms would 
then seem to be fixed by the field's corresponding fiat spacetime Lagrangians combined with 
the standard prescription for promoting them to curved spacetimes. 

As natural, and as straightforward, as the above prescription is for determining the cou- 
pling of matter to gravity, it nevertheless makes a number of implicit assumptions. When one 
tries to reconcile these assumptions with classical OR, a number of troubling inconsistencies 
become immediately apparent. 

The first implicit assumption is that the measuring apparatus does not play a role in the 
theory. That is, when calculating the various effects caused by the interaction of gravity with 
matter — such as, say, the scattering cross section of GWs — one does not have to explicitly 
include the measuring apparatus. This assumption is certainly true for all the other forces of 
nature; the existence of opposite or canceling "charges" for the EM, weak, and strong forces 
ensures that one can, in principle, screen out these forces, thereby making the measuring 
apparatus unaffected by them. This is not true for gravity, however; one cannot screen 
out the gravitational force. Why, then, should one not include the measuring apparatus 
explicitly in the construction of QG or of QFCS? Indeed, we shall show here that one must 
include the effects of gravitation on the apparatus in order to obtain physically correct 
results. 

One may argue that the scattering processes considered in quantum gravity occurs at 
very short length scales — the Planck length — and the presence of any measuring apparatus 
will have a negligibly small effect. However, we would expect from the correspondence 
principle that classical gravity results could be obtained — in some limit — from the quantum 
theory; indeed, the construction of QG 

mm 

makes explicit use of the classical theory. 
And in classical GR it is well known that the inclusion of the measuring apparatus — along 
with the observer — is crucial to understanding the dynamics of certain time- varying general 
relativistic systems involving tidal forces. 

Consider, for example, an isolated observer and a classical test particle initially at rest 
some distance away from him. While both the observer and the particle do not move spatially 
with respect to one another, they are both physical objects that move along their respective 
geodesies. If a GW now passes through the system, tidal forces will of course shift the 
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position and the velocity of the particle. However, these same tidal forces will also shift the 
position and the velocity of the observer; the observer cannot be isolated from the effective 
tidal forces caused by the GW. Thus, the observer cannot measure the motion of the particle 
independently of his own motion; he can only measure the relative motion of the particle 
with respect to himself. For GWs in the long-wavelength limit, the particle appears to the 
observer to undergo geodesic deviation motion [Eq. (35.12) of [Tf], and not geodesic motion 
as one might first expect. Indeed, a simple derivation of the geodesic deviation equation 
in this limit is to take the geodesic equation for the observer, subtract from it the geodesic 
equation of the particle, and expand the result in the distance separating the two geodesies. 

Based on this simple example we would expect that any physically measurable response 
of matter to the scattering of GWs calculated by using either QG or QFCS should include 
the effect of the GW on the observer. By extension, we would expect that in order to be 
consistent with classical GR, the construction of QG or QFCS should explicitly include the 
observer and his measurement apparatus from the very beginning. 

It may also be argued that, as with the other forces, explicit inclusion of an observer would 
be formally correct, but not required; the lack of its inclusion would not materially affect 
any calculation. This argument would also be in conflict with general relativity, however. 

Consider once again the simple system described above. When the geodesic deviation 
equations of motion are solved, one finds that the observed tidal response of the test particle 
to the passage of the long-wavelength GW is proportional to the distance separating the ob- 
server from the test particle; as long as the long-wavelength approximation holds, the further 
away the particle is from the observer, the larger its response to the GW. This ubiquitous 
response of classical matter to the passage of a GW is exploited in various GW detectors 
such as the Weber bar and LIGO (Laser /nterferometry Gravitational- wave Observatory); 
the larger the detector, the larger its response to the passage of a GW. The characteristic 
size L of the detector does play a role in the response of the system to the GW. 

Suppose that either QG or QFCS is used to calculate the response of a Weber bar or LIGO 
to passage of a GW through the system. It would be natural to use the complex scalar field 
^ to describe the system. Using the standard approach outlined above, it is straightforward 
to see that to lowest order, the coupling between a GW propagating in Minkowski spacetime 
with \1/ is ~ h^i^d'^'^'^d'^'^ . Even in the long-wavelength limit, this interaction term does not 
explicitly depend on the size of the system. Moreover, it is difficult to see how such size 
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dependence can be generated by this term. 

The second imphcit assumption made in Q and Q,S,Q is that one can always find a global 
time axis — and thereby construct a global coordinate system — in the curved spacetime. This 
is certainly possible for fiat Minkowski spacetime, which is often used as the background 
spacetime. It is also possible for the asymptotically fiat manifold that DeWitt considers in 
2]. However, we know from classical GR that it is not possible to find a coordinate system 
with a global time axis in general. 

The Minkowski spacetime and the asymptotically fiat spacetimes — along with the various 
black hole spacetimes — are stationary spacetimes. In these spacetimes one can always choose 
a frame where the metric does not depend explicitly on the time coordinate. Consequently, 
one can always construct a global timelike Killing vector, which can be used by all observers 
in the spacetime as their time axis (except, perhaps, at the event horizon or at an essential 
singularity). This Killing vector can then be used to construct what DeWitt termed the 
"preferred frame", or as Hawking and Ellis termed it, a "special frame" that all observers 
in that spacetime can agree to use. 

Timelike Killing vectors — and global frames — do not exist in general, however. Impor- 
tantly, they do not exist in the presence of a GW. Instead, each observer must choose his own 
local proper time axis, and construct his own local proper coordinate system from it. Conse- 
quently, one can only measure the relative motion between observers. This is the underlying 
physical reason why, in the example given above, one ends up with the geodesic deviation 
equation of motion |2| in the long-wavelength limit for GWs, instead of the geodesic equation 
of motion. 

n 

Although the above point is made very elegantly in the beginning of Chap. 4 of [8] for 

classical GR, it is relevant on the quantum level as well. As pointed out in both Chap. 3 

of and Chap. 3 of Q], the construction of Fock spaces for quantum fields in a curved 

spacetime is frame dependent; different choices of coordinates result in unitarily inequivalent 

Hilbert spaces. Thus, only for such spacetimes as the stationary and De Sitter spacetimes 

(considered also by DeWitt in where there is a "preferred frame", will it be possible for 

all observers to agree on what constitutes a particle state. It does not exist in general (see 
fl 

p] for a discussion of the relevance of the concept of "particles" in general spacetimes). 

As dissimilar as the two above implicit assumptions may be appear to be on the surface, 
they are nonetheless intimately connected. The experimental measurement of any physical 
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quantity requires an operational choice of origin, and a local orthogonal (tetrad) coordinate 
system. As any measurement is done through a physical apparatus, this mathematical choice 
of coordinate systems is fixed on a real physical object. The inclusion of the observer in the 
theory is thus equivalent to a choice of local coordinates; a choice of local coordinates must 
be equivalent to the inclusion of the physical observer. 

Although we have pointed out in the above a number of inconsistencies between results 
from classical GR on the one hand, and QG and QFCS on the other, the goal of this paper is 
not to present a reformulation of either QG or QFCS; we leave that task to future research. 
We instead address the issues raised by the two assumptions in the above by focusing on 
the dynamics of a much simpler system: the nonrelativistic, classical test particle in weak 
gravity. As simple as this system may be, especially when compared to the counterexamples 
we have listed above, many of the issues that we have raised above appear here as well. 
Fundamentally, what is at issue here is the appropriate choice of coordinates; this is an 
inherent aspect of classical GR, and is not due to a subtlety in the quantum theory. 

An analysis based on the dynamics of classical test particles has the added advantage of 
having limiting cases that have either been experimentally verified, or are in the process of 
being verified. In one limit, the Eotvos experiment, the advancement of the perihelion of 
Mercury, the deflection of light by the Sun, and the gravitational redshift are all calculable 
within the usual dynamics of test particles in stationary spacetimes based on the geodesic 
equation. In the other limit, the response of Weber bars and LIGO to the passage of GWs 
is calculable within the dynamics of test particles based on the geodesic deviation equation 
[3]. The result of the analysis in this paper must agree with these two limits; this serves as 
a stringent test of the validity of the approach we have taken and the coordinate system we 
have constructed. 

In the literature most analyses of the dynamics of test particles in curved spacetime are 
done in the same vein as the construction of QG and QFT in curved spacetimes, and are a 
direct generalization of the usual techniques for deriving Hamiltonians from Lagrangians for 
particles in fiat spacetimes. One starts with the usual geodesic action for the test particle 
moving in an arbitrary curved spacetime with a given metric g^,^. Time reparametization 
invariance of the action is broken either by choosing an explicit time coordinate, or by 
introducing a mass-shell constraint (by hand or through a Lagrange multiplier). Choosing 
x'^ as the general coordinate, the canonical momentum is calculated from the Lagrangian. 
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The Hamiltonian Hsf (SF for "standard formalism") is then constructed from this J9^ and 
the Lagrangian in the usual way. An analysis similar to this was followed by DeWitt in his 
1957 paper j^, albeit in much more detail, and in 1966 he applied the nonrelativistic limit 
of Hsf for charged test particles to the analysis of the behavior of superconductors in the 
Earth's Lense-Thirring field [lol |. 

It would seem that all we would have to do is to take the nonrelativistic limit of Hsf- 
However, the form of Hsf is dramatically different from the Hamiltonian for test particle 
motion derived in based on the geodesic deviation equations of motion. As with the 
case of QG and QFCS the same troubling questions come to the fore at this point: Where 
is the observer? What are physical quantities such as the position x and velocity v of the 
particle measured with respect to? What frame has been implicitly chosen by this analysis? 
Is this frame physical? We know from the observer-test-particle example given above that 
these are not fatuous questions. Rather, they directly address the underlying physics. 

It is certainly true that in some specific cases — such as the presence of a weak GW in 
the system — one can treat the time-varying part of the metric as a perturbation of one of 
the known stationary metrics; this time-varying piece would then be reflected as a pertur- 
bation on the Killing vector. One could then use the usual coordinate system for these 
spacetimes — augmented by the inclusion of the observer and his coordinate system — and 
calculate Hsf in the usual way. Doing so will not elucidate the underlying physics, how- 
ever, and it is difficult to see how the geodesic deviation equations of motion arises in this 
approach. 

The approach we shall take instead in studying the dynamics of nonrelativistic test parti- 
cles in the linearized gravity limit will be to construct a general coordinate system that builds 
in the essential physics from the very beginning. Since relative measurements between the 
observer and the particle always make physical sense, they are used as the foundation of our 
construction; the special case of stationary metrics will naturally be included. Specifically, 
we follow the considerations of [12] and jl^ (see also jl^ in connection with the coordinate 
system used in the analysis of LIGO): Every physical particle travels along a worldline fcir) 
with tangent vector (which does not need to be a geodesic) in the spacetime manifold 
M. Every measurement of the physical properties of the test particle by an observer must 
be done using an experimental apparatus. The observer — along with his apparatus — must 
propagate along his own worldline /"^(t) with tangent vector u^. Consequently, every phys- 
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ical measurement of the particle is done relative to the motion of an observer. In particular, 
in measuring the position of the particle, one measures the distance separating rdr) and 
/"^(r); in measuring the 4- velocity of the particle one measures of the relative velocity of the 
particle with respect to the observer 

Implementation of the above considerations proceeds quite naturally. As the observer 
prepares to take measurements on the test particle, he first chooses a local orthonormal 
coordinate system. In curved spacetimes, this involves the construction of a local tetrad 
frame . Naturally, this coordinate system will be fixed, say, to the center of mass of his 
experimental apparatus, and will thus propagate in time along the worldline /^^(t) as well. 
The observer uses the coordinate time of the physical apparatus to measure time, which, 
because he will not be moving relative to the apparatus, is also his proper time. Thus 
the time axis of the coordinate system he has chosen will always lie tangentially to FuIt). 
The position — which can be of finite extent — of the test particle is measured with respect 
to an origin fixed on the apparatus, and is the shortest distance between this origin and 
the particle. However, because the apparatus travels along its worldline, the origin of the 
coordinate system will also travel along a worldline in M. Later, when the rate of change 
of the position of the particle is measured at two successive times, the relative 4-velocity of 
the particle with respect to the apparatus will naturally be obtained. Thus, the observer 
constructs his usual laboratory frame that extends across his experimental apparatus, but 
now incorporating the nontrivial local curvature of M. We call this frame the general 
laboratory frame (GLF). 

Local coordinate systems fixed to an observer have been constructed before. The Fermi 

n 

normal coordinates (FNC) were constructed in the 1920's by Fermi [15[, and the Fermi 
Walker coordinates (FWC) were constructed in While an observer can use either set of 
coordinates, both make assumptions and approximations that drastically limit their useful- 
ness. The FNC — a direct implementation of the equivalence principle — are constructed so 
that the Levi-Civita connection F^^, vanishes identically along the worldline of the observer; 
only when one moves off the worldline does the curvature dependent terms begin to appear 
For the FWC, the restrictions on F^^, are somewhat relaxed, but certain components of 



only 

s 



F"^ — such as F^- where a, b are spatial indices — still vanish along the worldline. Once again, 
when one moves off the worldline curvature terms appear in the form of the Riemann tensor 
and its derivatives. In both, one effectively makes a derivative expansion in the Riemann 
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curvature tensor 
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In both FNC and FWC systems, choices for the value of F^^, — a gauge choice — have been 
made, and in both systems such gauge choices are inconsistent with the usual transverse- 
traceless (TT) gauge for GWs. While it is possible to study the interaction of GWs with 



test particles in these coordinate systems (see 



231 for FNC and |2J for FWC), doing 



so is cumbersome. For example, it has only recently been established that the TT gauge 
for GWs is compatible with the FNC |25^, but only in the long wavelength limit; the two 
are incompatible when the wavelength becomes smaller than the size of the experimental 
apparatus. In our construction of the GLF, no such restrictions on F"^ are made within 
the linearized gravity approximation. Thus, when we consider the case of GWs interacting 
with nonrelativistic particles, the TT gauge — or any other gauge — can be directly taken. 
Moreover, we do not make any restrictions on how rapidly the Riemann curvature tensor 
varies, and therefore are not restricted to only the long- wavelength limit. This enables us to 
study the effects of arbitrary-wavelength GWs on the motion of nonrelativistic test particles 
in large systems. 

It is here in our study of test particle dynamics that we obtain our most surprising 
result: Even though the underlying GW is a spin-2 tensor field, in the weak gravity, slow 
velocity limit, the GW acts on the particle through a local velocity field Na- This velocity 
field — which is an integral of the Ricci rotation coefficients — couples to the test particle as 
though it was a vector potential for a spin-1 vector field (see also ,20] for the additional 
terms that the Ricci rotation coefficients introduce in fermionic condensed matter systems 
and their implications), and its origin is the tidal nature of the forces that the GW induces on 
the test particle. It has the same properties as a vector potential: Like the vector potential 
for the EM field Aa, Na is a transverse field satisfying the wave equation. It is a frame- 
dependent field with the local Galilean group as its gauge group. Effective "electric" and 
"magnetic" fields can be constructed from Na in the usual way, and they are solutions to a 
set of partial differential equations that have the same form as the Maxwell equations since 
they are directly related to the electric and magnetic parts of the Weyl tensor, and thus to 
components of the Riemann curvature tensor. The equations of motion for the nonrelativistic 
particle have the form of a Lorentz force with the mass of the particle playing the role of 
the charge. As required, these equations reduce to the usual geodesic deviation equations 
0] in the long- wavelength limit. 
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The rest of this paper is organized as follows. In Sec. II we construct explicitly the GLF 
and its coordinates using a tetrad frame fixed to the worldline of the observer. The velocity 
of a test particle in the GLF is derived in the nonrelativistic limit. In Sec. Ill we use these 
velocities to construct the action, and then the Hamiltonian for the test particle in the GLF. 
We show that for stationary M this Hamiltonian reduces to DeWitt's Hamiltonian, and for 
ong-wavelength TT GWs propagating in a flat background it reduces to the Hamiltonian 



ll\ derived from the geodesic equations of motion. In Sec. IV, we study the properties 
of the velocity field A^^ introduced in Sec. 3 for arbitrary GWs, and construct effective 
electric and magnetic fields from it. These fields are shown to obey equations that have 
the same form as Maxwell's equations, and they are used to derive the equations of motion 
for a test particle. An Aharonov-Bohm-type interference effect for quantum test particles 
that is shown to follow from the effective vector potential Na can be found in Sec. V along 
with other concluding remarks. In Appendix A we present a brief review of the tetrad and 
linearized gravity formalisms, while in Appendix B, we derive the nonintegrable phase factor 
exp{i{m/h) f N^dX^} for the gravitational Aharonov-Bohm-type interference effect. 



II. CONSTRUCTION OF THE GLF 



As usual, g^v{x) is the metric on the curved spacetime manifold M with a signature 
(— 1, 1, 1, 1). Greek indices run from to 3, and they denote the coordinates for a general 
coordinate system on M. We will, however, be working primarily in one specific tetrad 
frame, and we will use capital Roman letters running from to 3 for the spacetime indices 
in this frame. (A summary of well-known results for linearized gravity and tetrad frames is 
given in Appendix A.) We reserve lowercase Roman letters running from 1 to 3 for spatial 
indices in the tetrad frame, and careted lowercase Roman letters for spatial indices in the 
general coordinate frame. A worldline with a timelike tangent vector parameterized by 
r is denoted as Fuir); this worldline needs not be a geodesic. Spacelike geodesies, with 
tangent vectors and parameterized by its arclength s, are denoted by 7x('S), and null 
geodesies with tangent vectors vr'^ parameterized by its arclength a are denoted by 7,1(0"). 

The construction of the GLF for the observer — being a specific choice of general tetrad 
frames that is fixed onto the worldline of the observer — is fairly standard. It must, however, 
be done without knowing the specific form of the underlying metric of M. Indeed, the local 
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metric at any given time is determined by making local measurements. We are aided in 
this construction by three observations. First and foremost, we note that the observer does 
not need a coordinate system that is nonsingular over all of M; such a coordinate system is 
known not to exist in general. All that is needed is a coordinate system that is nonsingular 
within the region of M where the observer makes experimental measurements. Second, we 
are working in the linearized gravity limit. This assures that we do not have to concern 
ourselves with coordinate singularities, and we can take curvature effects as perturbations 
on the flat spacetime metric. Finally, we are primarily interested in the effect of linearized 
gravity on nonrelativistic test particles; in this limit, incorporation of causality effects in the 
construction of the GLF simplifies dramatically. 

Let us consider an observer O with worldline ru{T). To perform experimental measure- 
ments at some time tq, O constructs a local orthogonal coordinate system centered on his 
experimental apparatus by choosing a tetrad frame {o^A^i'^o)}, ^ set of orthogonal unit vec- 
tors such that r)AB — o^a'^ i'^o) 9nu\ru{To) o^b" i'^o), where tjab is the usual Minkowski metric 
and gij,u\ry,{To) = o^AfiiTo) rj^^ o^BviTo) is the metric for M at r„(To). We use the presubscript 
o (for "observer") on to emphasize that at this point the frame only exists on /^^(r). 
Unlike the general tetrad frame, we require that u'^{tq) = o^o'^i'^o)', the time axis of the 
frame at tq lines up with the worldline of the observer. As usual, tetrad indices are raised 
and lowered by tjab, and general coordinate indices are raised and lowered by g^^. 

For the coordinate system at subsequent times we have to transport o^a'^ along FuIt) in 
such a way that that o^o always points along u^. If ^'^(t) were a geodesic, we would only 
need to parallel transport e^'* along it. However, because we are interested in nongeodesic 
worldlines we must instead use Fermi- Walker transport, a generalization of parallel transport 
that subtracts the nongeodesic motion of /"m(t) from the transport of o^a'^- ■^'^^ vector 
v'^ and a tangent vector to some worldline F^, the Fermi- Walker transport of v'^ along 

r^is 




(1) 



where as usual parallel transport along is 



(2) 



dr 



dr 



and is the connection on M. 
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FIG. 1: Parallel transport of ej^^ off ^'^(t) and along a spacelike geodesic 7j(-s). 

By the Fermi- Walker transport of oe^^(T) along Fuir), we find at each time r > tq, 

£%:M.0. (3) 

OT 

Not surprisingly qCq ^(t) automatically undergoes Fermi- Walker transport. The spatial 
tetrads, on the other hand, do not, and are solutions of the linear partial differential equa- 
tions 

= ^=%^-U»^^Ve„"(r), (4) 



dr V dr ^ 

with the appropriate initial condition at r = tq. 

To extend {oe^'"!''")} off Fu^r) we once again use the Fermi-Walker transport of e^^, 
but now in directions perpendicular to Fuir). At any fixed time r, let 7^(s) be a spacelike 
geodesic such that u^x^l = and 7^(0) = fu^r) (see Fig. [T]). We include the superscript 
T on 7^(s) to denote the implicit dependence of 7^(s) on r. For geodesies, Fermi- Walker 
transport is equivalent to parallel transport, and e^''(r, s) are solutions of 

0^^^, (5) 

along with the initial condition e^'^(r, 0) = oe^^(r) for each time r. It is straightforward 
to show from Eq. © that e^'^(r, s) eB^xij^s) = rjAB- Consequently, we can now consider 
to be a vector field on M such that 

Qfiu = e^/i (^Au, Vab = (^a' ^Bp.- (6) 
12 



While the above defines a set of orthonormal vectors for the observer, we still have to 
construct explicit coordinates for this frame. As mentioned in the Introduction, observer- 
based coordinates have been constructed before [3], [l^, although they were not explicitly 
constructed in a tetrad frame. Although the FNC and FWC are constructed using a series of 
approximations that make them unsuitable for our purposes, a number of the fundamental 
concepts used in their construction nevertheless carry over to our construction. In partic- 
ular, Synge [l^ introduces the notion of the world function to construct both coordinate 
systems (see also jl^). This extended object is a scalar function that measures the length 
squared between two points connected by a geodesic on M that are separated by a finite 
distance. It serves as a two-point correlation function that measures the net effect of the 
differences in local curvature between the two points. Moreover, because the world function 
is a length — and thus a scalar invariant — it is expandable in terms of the Riemann curvature 
tensor and its derivatives, thereby avoiding many coordinate-dependent artifacts. 

Useful though the world function may be, with the tetrad frame constructed above we 
have a more direct method of constructing coordinates. (This method is similar to the 

n 

approach followed in [20] for the FNC.) Like Synge, our method makes use of an extended 
object between two points on M, and we explicitly introduce a test particle V with worldline 
r'c(r) that is close enough to Fuij) for its physical properties to be measured by O's exper- 
imental apparatus. We ask what the coordinates of this particle are in (9's frame. To be 
consistent with O's experimental measurements, we parameterize /"c(t) by the proper time 
of (9, not V. Then let X^(r) be the position of V at any time r in the observer's frame. 
X^{x) : can also be considered as a coordinate transformation from the general 

coordinates to the tetrad frame at any time r, which in a small neighborhood Uo of O 
is 



dx'^dx^, (7) 

Uo 



Uo 



SO that from Eq. (jH)), 



Taking the derivative of the first equation in Eq. (jH)) with respect to x'^, we find that the 



integrability condition: 0,^6^ = d^e^. This condition only holds within Uq (see Eq. ()A3 
To extend X^ off O's worldline, we note the following. 
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ul'(T) 



FIG. 2: A sketch of the worldlines /^^(t) and rc{T) of observer O and test particle V, respectively, 
with the spacelike 7^(s) and null 77r(o") geodesies used to construct X"^ shown. The end points of 
both geodesies are fixed onto specific points on Fuir) and rc{T). The simultaneity with respect to 
the observer O in time along 7j(s) is an approximation that only holds in the linearized gravity 
and nonrelativistic limits. 

Solutions of Eqs. ^ are clearly path dependent. For spatial components X^, we consider 
again the spacelike geodesic 7^(s), but now connecting -r„(r) to Fc^t) such that 7^(0) = 
r„(r) and 7j(si) = rd^r). (See Fig. El) By integrating Eq. (jHl) along this geodesic, we obtain 
the spatial coordinates of the test particle 

/•7j(^i) f 



as a straightforward extension of the tetrad framework to coordinates. We have made use 
of differential forms through ba = CA^dx^ (see Appendix A); this will greatly simplify our 
analysis later. In defining Eq. Q, we have explicitly assumed that the partial differential 
equations in Eq. (jHJ are integrable. While this is not true for general spacetimes because 
of the presence of singularities, they are integrable in the weak gravity limit that we are 
working in. We also note that the length of Xa is proportional to the proper length si of 7^. 

Like the world function, Xa(r) is an extended function of two points — one on the worldline 
of O and the other on the worldline of V. Indeed, it is straightforward to see that the spatial 
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coordinates in the FWC are an approximation of Eq. (jH)) by taking ^is a constant equal 
to its value on 7^(0); the remaining integral is proportional to the world function. We 
emphasize, however, (the index a notwithstanding) that Xa is the integral of a differential 
form, and is a scalar function on M (see Appendix A). 

The construction of the time component of the test particle Xq{t) is more complicated 
because of causality. As Synge pointed out in his derivation of the FWC, using space- 
like geodesies 7^(s) in Eq. Q is somewhat artificial; no physical measurements ever take 
place along spacelike geodesies. Strictly speaking, we should have instead used null vec- 
tors — corresponding to optical measurements — in the above. However, this would have 
resulted in a set of optical coordinates, and since we are primarily interested in the motion 
of a nonrelativistic test particle, would have been needlessly complicated. Instead, we note 
that in the nonrelativistic limit the forward lightcone of the observer opens up, and the 
observer's null geodesic is well approximated by a spacelike geodesic in this limit; we would 
thus expect the above construction of Xa to be valid in the nonrelativistic limit. The same 
argument cannot be made for the coordinate time Xq of P, however; causality still has to 
be taken into account. For this coordinate, we first construct Xq using null geodesies, and 
then take the appropriate nonrelativistic limit. 

Figur^ shows explicitly the null geodesic and spatial geodesies that we use in this con- 
struction. At a time r' < r, let 7,1(0") be a null geodesic that connects O at time r' to V at 
time r: 77r(0) = fu^r') and 7,1 (o"i) = Fc^t). We define 

Xo{r) = '\oJr{cr),s{a)]7i^{a)da- (- [ eo^(f , 0)u^(f )rff ) . (10) 



^7vr(0) V / \ Jt' / 

The first term in Eq. (fTUI) is the time it takes an optical signal to reach V. The second term 
is the amount of time that passes for the observer for the optical signal to reach V. Since 

Xo(r) = -{r-T')+ [ eo. (11) 

Note that because u^Xn = 0? unlike Xa we cannot simply replace tt^ by the nonrela- 

tivistic limit; the second term in Eq. (jllj) would vanish automatically. This limit has to be 
taken much more carefully 

The position of the test particle was arbitrary and could have been placed at any point 
near /^^(r). Thus the GLF is a combination of a tetrad frame {e^'^} fixed to the worldline of 
the observer together with the coordinates (r, X"). It is important to note that X'^ measures 
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the relative separation between the worldhnes of O and of V. For "small" X°' 39|, and for 
r„(r) and Fc^r) geodesies, X'^ is simply the geodesic deviation between O and V. We have 
also chosen, as most physical, to use the proper time of the observer as our time coordinate; 
Eq. (jlip gives the dependence of the test particle's coordinate time on r. As we usually 
use the frame where the observer is at rest in the GLF, r coincides with the coordinate 
time of O. We shall simply use r as the time variable for O to avoid introducing additional 
notation. 

In what follows, the final expressions of all physical quantities measured in the GLF 
will be expressed in terms of (r, X"), the coordinates that O measures in the GLF. In the 
frame, = C^^m^ = 6q, and u"^ points along the time direction for the observer, while 
X^i^, s) = e"^ ^x^{t, s) is a scalar function in the GLF such that 

The first term vanishes because 7^(s) is a geodesic, and the second term vanishes from the 
construction of e^ ^ in Eq. 0. Thus in the GLF = X^i"^) only. In addition, since 
= Xm^'* — Xa^q^ X° = and x^ is a unit spatial vector pointing directly at the test 
particle at any time r j4Q']. As a final point on notation, although da = d/dX"- = e^^^d^, 
to avoid confusion we shall always write d/dr instead of do, which is instead reserved for 
d/dX^. 

We now turn our attention to finding the 4-velocity of a test particle as measured in the 
GLF. To do so we refer to Fig. IHl which shows the position of both the observer and the 
particle at two subsequent times r and t + 6t. Note that both the observer and the observed 
are moving: O along its worldline /^^(t) and V along its worldline Fclr). The observer can 
only measure the relative 4-velocity between V and himself. Beginning with the spatial 
coordinates, we take 

dX, _ Xa{T + 6t) - Xa{r) 



dr St^o St 

ea}, (13) 




where 7J!+'^'^(s) and "yZ{s) are spacelike geodesies from O to V at t + 5t and r, respectively 
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FIG. 3: The location of O and V at two subsequent times, and the spacehke geodesies used in 
constructing Xa as they propagate along their worldlines. The closed rectangular region D used 
in Eq. (|14|) is shown also. Notice that the vector ^'^ as it varies along 7j(s). 

(see Fig. 3). Adding and subtracting integrals along the worldlines of O and V, 



dXg _ 1 
dr 



Uni ^ W + y 0)u^\'r)df - J ea/,(f, Si)c^(f)cif j , 

Ca{T,X-)- lim ^ f dea, (14) 



St^o St Jd 

where ca = eA^i^T, X'^)c'^{t) is the 4- velocity of the test particle in the GLF. D is the closed 
region bounded by 7^(s), 'y^.'^^'^ (s) , and the worldlines /^^(t) and Fc^r) between r and t + 6t. 
We parameterize the region D by the 1-forms d/^ = x^ds and d/2 = ^'^dr where 

?" ^ ^. (15) 

is a function of (t, s). Because x^X'* = I5 then ^^x'^ = 0^ ^ind d/J* ^'^d ^^^2 linearly 
independent. From Eq. ()A3|1 . de^ = A UaB where cja^ = e^'^V ^esudx^ is the Ricci 
coefficient 1-form. Then, 

/ de„ = - / / {e'^^ujjyaB - e^^cu^afi) ^^x^'dfds. (16) 
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The limit 5r ^ is now trivial to take, leaving only a path integral along 7^(s). Then 

Ca = —r^+ / {uJBac- (^caB}^^dX'', 

dr Jo 

= ^ + l {(e,^ac)e^-(e,eVaB}, (17) 

where we have now expressed the path integral in GLF coordinates by using dX^ = ^X^df, 
and = e^^^^. As usual, ujac) = ^^^^f^ac- 

To determine ^, we note that measures the deviation in the geodesic 7^^''^(s) from 
"iZi^s) at any time r along s, and is thus the solution of the geodesic deviation (or Jacobi) 



equation in s 



121. In the GLF, 



= ^ + R\e (18) 

where = and ^ = cbdX^X^ (see Appendix A). Note that x^ is independent 

of s and only fixes a direction in the above. 

From Fig. El we see that interpolates between the tangent vector to the worldline of 
C, and the tangent vector to the worldline of V. Since ^"^Xa = i^Xn = 0, the precise 
boundary conditions are: lims^o+ 'C^('i") -5) = and lim^^^- ^^(r, s) = {S^ — x^Xb)c^- 
However, appears in Eq. through of the 2-form d/^ A dl2. Since dZg li^s parallel to 
X^, this automatically projects to zero any component of parallel to x"^. Consequently, 
we can without loss of generality replace the second boundary condition by the much simpler 
condition lim^^^- ^^{t, s) = c^. Equation (fTH|) can then be solved iteratively, and 

^ir, s) = - [sc^ir) + {s, - s)6^} + H ^{r, s')e{T, s')G{s, s')ds', (19) 
■51 Jo 

where G{s, s') is the Green's function 



s , , . . s' 



G{s, s') = — (si - s')9{s' - s) + — (si - s)e{s - s'), (20) 

Si Si 

and 6{x) is the Heaviside function. 

For Co we proceed the same way. Using now the diagram in Fig. HI 

Co = ^ + lim ^ [ e'^A u^oa- (21) 
dr St-^o ot ' 



D' 



To parameterize D', we again take d/f = ir^da, but now dli^ = C^dr where, in the GLF, 

(^{T,a) = 6^ + {c^{T)-6^}^^. (22) 
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z') 



FIG. 4: The null geodesies used in constructing Xq at two subsequent times as O and V propagate 
along their worldlines. The closed region D' used in Eq. (|21j) is shown. 

Clearly, ("^ and are linearly independent. Unlike the spatial components of the 4- velocity, 
we cannot use ^Dn^^ / dr to parameterize D', because it is either a null vector or a spacelike 
vector. Although like C,^, interpolates between the tangent vector at Puir') and the 
tangent vector at ^'^(t), it does not include the corrections due to nonvanishing curvature 
that Eq. p8|) does. Because we will be working in the linearized gravity limit, such correction 
terms are of the order Hab times the particle velocity, and can be neglected to lowest order. 
We therefore get 



Equations ()17|) and (j^Hj) in principle determine the components of the 4-velocity of 
the test particle in the GLF. However, because both and themselves depend c^, these 
equations form a set of coupled integral equations in c^. While these equations can be 
solved iteratively using Eq. and Eq. ()22|1 . we are primarily interested in the behavior of 
nonrelativistic test particles in the linearized gravity limit. In fact, much of our construction 
of the coordinates for the GLF is only valid for the nonrelativistic test particle. We thus 
keep only terms linear in h^^ where g^i, = rj^^ + h^^ (see Appendix A), and we approximate 
c'' ~ 1, keeping only terms linear in c"^ — 5^. Since both the spatial velocity and the 
curvature effects are small, we also neglect cross terms of 0((c'^ — 5^)/iA_B)- Thus, we can 




(23) 
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neglect curvature corrections in Eq. ()19|) altogether, and can take ^ Sq since always 
appears in the combination C,uj. Using a similar argument, we take ~ Sq as well [which 
is why we did not have to be concerned with curvature corrections to ("^ in Eq. (|22p]. With 
these approximations, 

dr 



dX 

CA = -T^+ I WoAb - UJbAo} dX\ (24) 







where dX^ = x'^dr. 



III. THE HAMILTONIAN 



The action for a test particle with mass m and charge q in general coordinates is 



I = -m J ^J-Cf,c^^dT + q j A^c^dr. (25) 

Because Eq. is time reparameterization invariant, we are free to choose the proper time 
r of (9 as our parameterization. In the GLF the Lagrangian becomes 



£ = -m^/id^y - (c«)2 + qAqC^ + qAaC", 

1 /"dX" Z"^" ~ \ ^ f^" 

^ -m + -mf-^+ / H%-cU;,%}rfxM -m/ uJ^^^dX^ 

( Z"^" \ ( dX"" Z"^" ~ \ 

+gAo f 1 + y ^oobrfxM + qAa [—^ + - rfxM , (26) 

where Ab = e^^A^ is the vector potential in the GLF. In the above we have used the non- 
relativistic and linearized gravity limits. This includes taking dX^/dr ~ 1. The momentum 
canonical to X"" is then: 

5C dX f^'^ 

5{dX^/dr) ='^^^'^L {-o..-^.ao}rfX^ + gA., (27) 

so that in general 

-qAo(^ + tuoobdX'^ + m uooohdXK (28) 
There are two special cases to consider. 
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A. Stationary metrics 



For stationary metrics there is always a frame where gf^i, is independent of time. In this 
frame, 

^OOb = -dbhoo, Uoab — i^baO = -dbhao, (29) 



where we have used Eq. ()A9|) . Then 



d X 1 

Co = — ^ + -{/ioo(r,X'^)-/ioo(r,0)}, 
dr 2 

dX 1 

Ca = ^ + -{/iao(r,X")-/^,o(r,0)}, (30) 



while 



H= ^{Pa- qAa)' - -Pa {Ki^, - K{r, 0)} + -m {/ioo(r, X") - /ioo(r, 0)} 

-qAo {1 + /ioo(r, X") - /ioo(r, 0)} . (31) 

To compare this with DeWitt's results 0], we first remember that we are working in the 
(orthogonal) GLF, while DeWitt is working in the general coordinate frame. To transform 
the above back to the general coordinate frame, we use Eq. ()A8jl in 

(1 \ Tn 
sil-^Kjp,^Pa-jh^, (32) 

where we follow DeWitt and neglect terms of 0{ph). For clarity, we use a caret to distinguish 

between indices in the GLF and DeWitt's frame. Like DeWitt we neglect terms 0{Ah) as 

well, and find 

1 / in \ ^ in 

H=^[Pa- mK^{r,X^) - -/i,o(r,0) - qA,) +- {^oir^X^ - /ioo(^> 0)}+gAg. (33) 

This agrees with DeWitt's result up to a constant shift in velocity and in energy. This shift 
is needed because DeWitt's coordinates are fixed to the origin of the mass generating the 
gravitational field (for example, on the center of the Earth), while our origin is fixed on the 
observer (on the surface of the Earth, say). 



B. For gravitational waves 



For GWs we work in the usual TT gauge, then cjoob = now, while 
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from which we see that 



^0^^, c. ^ - iV„(A-,, (35) 



where 



iV„(X") = / UJaO. (36) 

Then, 







H = 7^{Pa- qAaf + PaiV" " qA^. (37) 
2m 

Now, in the long-wavelength limit, Na = X^^dhab/dr) /2. Thus, for neutral particles Eq. ()37|1 
agrees with the Hamiltonian derived in from the geodesic deviation equation. In that 
paper one of us (A.D.S.) also postulated the existence of a minimal coupling with the EM 
field which was of the form ~ {pa + Na — qAa^ using the notation of this paper. This was 
based on the usual arguments for minimal coupling, and was valid within the framework 
and approximations of that paper. This direct coupling between Na and Aa does not appear 
in this more general analysis. However, a direct coupling betwen the spin connection and 
the EM vector potential is present for fermionic systems as presented in j26,j . 



IV. PROPERTIES OF A^^ 

If we expand out the kinetic part of the Hamiltonian Eq. ()37|) . we see that for small Aa 
both Aa and A^^ couple to the test particle in much the same way. Surprisingly, this simi- 
larity between Na and Aa goes much deeper. Like Aa, Na — which has units of velocity — is 
transverse and satisfies the wave equation in GLF. Like Aa, Na is a gauge dependent object, 
the gauge symmetry in this case being the choice of local Lorentz frames, and the gauge 
group being the Galilean group (since we are dealing with nonrelativistic test particles). 
As with Na, we can construct from Na effective "electric" and "magnetic" fields that obey 
partial differential equations that have the same form as Maxwell's equations. These fields 
have direct physical meaning, and they turn out to be proportional to the integrals of the 
electric and magnetic parts of the Weyl tensor. 

Consider first 

dN 1 

— ^ = \im-{Na{T + 6T,X'')-Na{T,X'')}, 

or St^o ot 

= cuB,„(r, X")c^(r) -cuB,o(r,0)5o^- lim ^ [ dcjao, (38) 
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FIG. 5: A sketch of the two spacehke geodesies 7^(s) and 7^/(s) for X" and X"' -\-6X^ respectively 
used in calculating spatial derivatives. The addition of 7j^(s) linking the two points defines the 
closed triangular region A used in Ea.(|^T|). 

where we have used the same arguments as in Sec. II. Once again, because we are in the 
nonrelativistic limit, the domain D' in Fig.|3]is well approximated by the domain D in Fig. 01 
Using Eq. ()A4|) and the usual parameterization of D, 

lim — / RBCaodlf A dl^, 
D 

RbOaodX''. (39) 

Jo 

since ~ 5q and c^oao = for GWs in the TT gauge. Similarly, 



dr 2 St-^o 5t 



d'^N 1 C 

lim — dcRsoaodl^ A d/| 



B 

5 



doRboaodX\ (40) 







since i?oo6o = 0. 

For the spatial derivatives, we refer to Fig. El and now consider the three spacelike 
geodesies bounding the closed surface A formed by 7^(5), connecting the origin to X"-; 
7^,(s), connecting the origin to X'* + 6X^ for 6X^ small; and 7j^(s), connecting X'* to 
X'' + 5X\ Then, 
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SX^0 6X~b ' ' '^aO- / d^aO ) , (41) 



A 



where we have once again added and subtracted an integral — along 7^^(s) — and used 
Stokes's theorem on the closed boundary dA of A. Parameterizing A by dl'( = x^ds 
and = x'^ds', then 

dbNa = iOtaO - / RcbaodX'. (42) 



^0 

This equation shows explicitly the difference between taking spatial derivatives in flat space- 
time verses taking derivatives in curved space. In flat spacetimes the gradient of A^^^ would 
simply be Ubao- In curved space, on the other hand, we have to take into account the differ- 
ences in curvature between X"" and X"" + 6X'', and that introduces the additional curvature 
term. Proceeding similarly, 

dbdcNa = dcUJbao - Rcbaoir, X") + ^hm^-^ j deRdbaodl^ A dl^, 

= dcUJbaO- Rcbao{T,0) - / dcRebaodX" . (43) 

^0 

From Eq. (021), 

d'^Na = u;\, + / RaodX''. (44) 



For GWs in the TT gauge cu'^^q = ((9/i^/9r)/2 = and Rao = -□/iao/2 = 0, where □ is the 
d'Alembertian operator. Na is thus transverse, and d"'Na = 0. Next, from Eqs. and 

= - (^) (45) 







Or , 

But this also vanishes since Rab = 0. Thus, like the vector potential for EM, A^^^ is a 
transverse vector field obeying the wave equation □A'^q = 0. 

Next, to construct the effective "electric" and "magnetic" fields from Na-, we consider 
the four- vector A^^, constructed from Na with A^q = 0, much like the vector potential for 
EM in the Coulomb gauge. Then clearly O^Na = 0, and Na is a transverse field. We then 
define an effective field strength ^ab = QaNb — ObNa with the corresponding "electric" and 
"magnetic" fields 

dNa 



aO 



Or 



^ l^a6c^^^ ^ e^^'^S.AT,. (46) 
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Using the same methods as above, we see that 



Or L V dr 



hOaO 



dX" 



db<Ea = RbOao{T,Q)+ / dbRcQaodX", 

Jo 



e 



dr L V dr 



d,^- = -e-f-R,f^^{r,0) - [ e^f'dtRefcodX'', (47) 

Now, it is clear that by definition e'^^'-^^dsdcD = 0, and this relation must explicitly 
hold if ^AB has been defined consistently. To verify this, we first consider 

d'^^a = -e'^f'Rafcoir, 0) - / e'^f'daRefcodX' = 0, (48) 

Jo 

which vanishes identically using the 1st (Eq. ()A5jl ) and 2nd Bianchi Identities (Eq. ()A6j) ). 
Similarly, 

^"^^ - e^'^d^e, = e'^'^ r dX^ (_^^ + d^R^ , (49) 



or Jo \ ^'^ 

also vanishes for the same reason. Thus, e^^''^^ Ob'^cd = is as a direct consequence of the 
Bianchi Identities for the Riemann curvature tensor. 

For the other part of Maxwell's equations, d^^AB = 0, we consider 

d^(Ea = -Rooir, 0) + / d'^R.oaodX'^. (50) 

From Eq. (jX6|) . 

= d^RcooA + d.Roo - (51) 
Once again, Rab = implies that = 0. Finally, we consider 

e-a<B. + ^ = -r7^VA- = 0, (52) 



dr Jq \ dr 

using Eq. ()A6|) and -Rab = once again. Consequently, (£a and 05" satisfy equations that 
have the same form as Maxwell's equations. Notice also that like EM waves, d^^AB = 
only in the absence of sources. 

The surprising connection between (Ba and OS", and Maxwell's equations can be under- 
stood by looking at the "electric" and "magnetic" parts of the Weyl tensor Cabcd- For 
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GWs propagating in a flat background, Cabcd = Rabcd and the Riemann curvature tensor 
can be separated into two parts The electric part, 

Eab = —CoaOb = —RoaOb, (53) 



HI = ^e^^'fCefob = -^e'^'^Refw. (54) 



b ma _ / ijajvb 



and the magnetic part. 
Clearly 

= - / EabdX", ^- = - I H^dX". (55) 
Jo Jo 

(£a and 25" are simply path integrals of Eab and ; they obey Maxwell's equations because 
Eab and obey tensor Maxwell-like equations |l3|. 

If Na functions as an effective vector potential for the GW, what, then, is the correspond- 
ing gauge group? Notice that although we have chosen a speciflc coordinate system — and 
thus broken general coordinate invariance — we still have a residual invariance left over. To 
see this, let us do a local Lorentz transformation L^"^ on e^^ such that tiab = Lj^^L^^ 
while e^'^ = L^^e^^. This leaves g^j^i, = ca^i y = e."^ ^ invariant, and we see that local 
Lorentz invariance is still left over. Now, it is straightforward to show that under a local 
Lorentz transformation ujcab ^cab + ^^^c^da'j ^cab transforms anomalously under a 
local Lorentz transformation. Consequently, 

Na^Na+ / L^'MoadX^ (56) 







Because we are working in the nonrelativistic limit, the gauge group for Na is the local 
Galilean group. Indeed, for a pure boost in the nonrelativistic limit, L^^ ^ — VaK°^^ 
where -f^^^ is the generator of boosts. Then, 

Na -^Na- dv-a = iV, - {t;,(r, X'') - Vair, 0)} , (57) 

^0 

and Na changes by a local velocity fleld. 

Finally, we calculate the equations of motion for the test particle in a GW. Equation 
for GWs is 

from which we get 

™^ = -™ (S" + + 4^" + ^""^ftl . (59) 



dr'^ 1 dr \ \ dr 

26 



after using Eq. ()46p : like DeWitt we drop terms 0{Ah). In Section II we stated that we 
would neglect terms of order (c^ — 6Q)hab, and yet Eq. contains just such a term. It 
is, however, straightforward, though tedious, to repeat the calculation for Eq. keeping 
the next higher order velocity-connection terms. After doing so we still obtain the above to 
lowest order. 

Let us suppose that a GW can be approximated as a plane wave with wave vector ka- 
Notice that the integral in Eq. is always over the plane perpendicular to ka, yet hab 
is a only function of the spatial coordinates parallel to ka- Consequently, Na = UbaoX^, 

= RbOaoX^ and 03" = e^'^^RbceoX'' /2. Thus, for planar GWs, 

= [Rma,X' + '^X^Rabco^ +q\Ea + ^abc^B^^ . (60) 

As expected, in the low velocity limit where the 03" term can be neglected, this is just the 
geodesic deviation equation for a charged test particle. Note also that Eq. (jUUj) holds for 
general GWs in the long-wavelength limit as well, and that it agrees with the equation of 
motion for a charged test particle interacting with a GW found in 



V. CONCLUDING REMARKS 

In Eqs. fl26p and ()28|) we can clearly see the natural separation that occurs in the dynamics 
of test particles in two extremes. In the one extreme the metric for M can be approximated 
as being stationary. Matter moving with a characteristic velocity v much smaller than the 
speed of light {v/c << 1) makes up the dominant contribution to the stress-energy tensor in 
these spacetimes; thus to a very good approximation the presence of GWs can be neglected. 
In this extreme our Hamiltonian reduces to DeWitt 's Hamiltonian, and the straightforward 
approach to the derivation of Hsf described in the Introduction works. In the other extreme, 
all the curvature effects are due to GWs, and the characteristic velocity for the metric is the 
speed of light. In this extreme, our Hamiltonian reduces to the Hamiltonian derived in 
for geodesic deviation motion in the presence of GW's in the long- wavelength limit, and the 
equations of motion reduce to the usual geodesic deviation equations of j7|. 

Let us be very clear. The dynamics of nonrelativistic, classical test particles in stationary 
spacetimes that we have derived in the GLF and are the same as the dynamics for these 
particles derived using standard methods. Thus, such classical tests of GR as the perihelion 
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of Mercury and the gravitational redshift will follow through in the GLF as well. By ex- 
tension, because QFCS's are usually formulated in stationary spacetimes, we would expect 
such an analysis as the evaporation of black holes and black hole thermodynamics to hold 
also. It is when we consider nonstationary spacetimes — such as those where GWs play a 
dominant role in determining the physics — that we must take care to explicitly include the 
observer and his experimental apparatus in the analysis. 

Traditionally, there is an almost linear progression of approximations, such as the hi- 
erarchy sketched in Fig. IHK) that we usually associated with GR. Within full GR we can 
make the linearized gravity approximation, which is thought to be inclusive of both the 
)arametrized-post-Newtonian [7| (stationary) and the GW (nonstationary) limits (see also 
). We can, with further restrictions, then pass over to either the parametrized- 
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post-Newtonian (PPN) limit or the GW limit within this linearized approximation. Based 
on the arguments above and the results of our analysis in this paper, we find the separation 
between the PPN and GW limits to occur at a much higher level; the dynamics of particles 
in stationary spacetimes is drastically different from that in nonstationary spacetimes. This 
conclusion is consistent with 0, and 0, all of whom note differences between dynamical 
theories in stationary versus nonstationary spacetimes. Thus, instead of the standard formu- 
lation consisting of the linear sequence of approximations in Fig. we should, as shown in 
Fig. inb; instead separate the dynamics in stationary versus nonstationary spacetimes from 
the start. 

What is most surprising about our results is the form that the coupling of general GWs 
to the test particle takes in the GLF. The tidal nature of the force that a GW induces on a 
test particle introduces an effective velocity field Na to the system that couples vectorially, 
under the Lorentz Group, to the particle, even though the GW itself is a rank-2 tensor. 
This Na — being a path integral — is in general a nonlocal function of the of Ricci coefficients, 
although for planar GWs it reduces to the product of the Ricci coefficient and the position 
of the test particle. [A vector coupling of GWs to matter was proposed by one of the present 
authors (R.Y.C.) in |29] based on a generalization of the PPN formalism.] Thus the motion 
of the test particle is determined by the cumulative effects of the GW along the distance 
between the observer's worldline and the test particle's. Like the vector potential for EM, 
Na is a transverse field satisfying the wave equation, and like EM, effective electric € and 
magnetic ?B fields can also be constructed from Na that obey equations that have the same 
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FIG. 6: Fig. 6a is a schematic of the hierarchy of approximations that we usually associate with 
GR. Fig. 6b is the modified hierarchy based on the results of this paper. Notice the explicit 
separation between stationary and nonstationary spacetimes from the start. 

form as Maxwell's equations. Thus, Na functions as a vector potential for € and 03 with its 
gauge group being the local Galilean group. Indeed, the equations of motion for the test 
particle in the presence of arbitrary GWs has the form of the Lorentz force with the mass 
m of the particle playing the role of the "charge." 

Surprising though this vector coupling may be at first glance, after further study we see 
that there are deep connections between A^^ and the tetrad formulation of eneral relativity. (B 
and ^ are proportional to the path integrals of the electric Eab and magnetic Hab parts of the 
Weyl tensor, respectively, and it is precisely the Weyl tensor that contains GW excitations. 
It is also known that Eab and Hab obey tensor Maxwell's equations for GWs propagating 
on a flat, source-free background, and the fact that € and 05 satisfy Maxwell's equations 
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is simply the reflection that Eat and Hat do. In fact, the identity e'^^^^dsdcD = holds 
precisely because of the 1st and 2nd Bianchi conditions for the Riemann curvature tensor, 
while d'^'^AB = only holds because we are in a source-free region for GWs. As for A^^^ 
itself, it is precisely "half" of the loop integral of the real part of Ashtekar's connection in 
the loop- variable formulation of quantum gravity. Indeed, Ashtekar's loop integrals can be 
thought of as quantum corrections to the "classical" Na presented here j30|. 

Since Na is a nonlocal object, it is natural to question the practicality of Eq. as 
compared to, say, Eq. ()60|). whose form for long- wavelength GWs is well known. Indeed, 
it would seem that for Earth-based systems the long-wavelength version of Eq. ()fi()|l would 
be sufficient to describe all physically interesting systems. We note, however, that at 4 
km the arms of LIGO are comparable to the reduced wavelength of the GWs for GWs 
with a frequency of 10 khz, which is at the upper limits of LIGO's frequency response 



spectrum [31[; measurements of the spatial variation of GWs are becoming experimentally 
accessible. Similarly, when LISA's (Laser Interferometer S'pace Antenna) arm lengths begins 
to approach the reduced wavelengths of GWs in the high-frequency limit of its sensitivity, the 
concept of the Na field and the GLF becomes essential. For these arbitrary GWs, Eq. 
instead of Eq. (jUUI) must be used. Since it is Na in general that will be measured, and not 
hab, it is an interesting, and important, question as to how much information about hab can 
be obtained from measurements of A^^. 

While most of our focus has been on classical dynamics, we end this paper by looking 
at a gedanken experiment that probes the essential difference between the classical and 
quantum dynamics of nonrelativistic test particles in the linearized gravity limit. Extension 
of Eq. (|?fj) to quantum mechanics is straightforward. For nonrelativistic, neutral particles 
with wave function ip, 

th^ = -I^V^V; - ihNad''^. (61) 

Notice that the interaction term Nad""ip has the same form as the coupling of an EM vector 
potential Ab (expressed in GLF coordinates) to a charge particle coupled to weak EM 
fields. In the EM case, the solution of equations of this form leads, on the quantum level, to 
the Aharonov-Bohm effect, which is expressible through Yang's nonintegrable phase factor 
(Wilson loop) 2^ Xem = exp{i{e/h) f ABdX^}. In direct analogy, Eq. (jF)T|l suggests that 
like the EM case, a time- dependent Aharanov-Bohm interference experiment such as those 
described in jsSi or js^ can in principle be done, but now with GWs. We would expect this 
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to lead to a corresponding phase factor for gravity Xgr = exp{i{m/h) ^ NAdX"^}. While 
actually performing this experiment would be unrealistic at this time, the corresponding 
gedanken experiment does illustrate both the essential difference between the classical and 
quantum systems, and the importance of this gravitational nonintegrable phase factor. 

Consider the gedanken Aharanov-Bohm-type experiment shown in Fig.[7|using a Gaussian 
GW packet with a width low at some time tj. Fig. [7^ shows the spacetime diagram for 
the particle and the GW used in the interferometry experiment with a time slice drawn in 
at t/. A schematic of the corresponding physical apparatus at this time slice is shown in 
FiglZb- We require that the GW also be an unipolar wave packet so that at all times the 
amplitude of the <E and 53 for the packet is non-negative. The interfering particle is taken to 
be a Gaussian wave packet as well, but with width Ip at any time t. After being emitted at 
the source, the particle at event pA passes through the first beam splitter shown in Fig. [7|d, 
and there is a finite probability amplitude it will propagate either along or along F"^ 
until it is is recombined at event pb- It is impossible, even in principle, to know which path 
the particle took through the interferometer. As usual, the combined path F = F^ U F"^ 
encircles the GW beam. 

The Bonse-Hart interferometer of a type shown in Fig. [7|d is a concrete example of an 
interferometer that could be used in this experiment; in this case the interfering particle 
could be neutrons. The gray oval patch at the center of the interferometer represents the 
time slice at tj shown in Fig. [7^ of the Gaussian beam waist. Most importantly, we choose 
the size of the interferometer such that R » lew + where R is the distance from the 
center of the GW beam in Fig. UJp to each arm of the interferometer. Thus, on a classical 
level the GW and the neutron worldlines do not intersect, and therefore the neutron feels 
no forces at any time arising from the GW. On a quantum mechanical level, the amplitude 
of the neutron's wave function is exponentially small where the amplitude of the GW is 
large, and the amplitude of the GW is exponentially small where the amplitude of the wave 
function is large. 

Although the wave function of the neutrons satisfy Eq. (pTT|) . because they are Gaussian 
wave packets we can make a WKB-like approximation by taking ip ^ e*®. For hat with small 
spatial variations along the worldlines, Eq. reduces to 

= ^ + Nad^Q, (62) 
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where we have neglected terms (9(|VGp) since Na is small. The solution of Eq. ()62p in the 
linearized gravity limit for the particle propagating along worldline is 



^r^ = J I NAdX^, (63) 



from Eq. ()B5|) of Appendix B. Similarly, if the particle had traveled the worldline F 



2 



6^2 = ^ / NAdX"^. (64) 

r2 



h 

The phase factor 



r Tn f 

exp{iAe} = explGj^i - Qp^} = exp J NAdX' 



(65) 



is a measure of the phase difference between the particle propagating along F^ versus prop- 
agating along F"^. Equation is simply the \gr that we predicted above. Next, since F 
is a closed loop, by the Stokes's theorem, 

exp{iAe} = exp j ^AsdX^ A dX^j , (66) 

where V is the surface bound by F and we have used the definition of ^ab- Because the 
boundary of V is made up of timehke curves, from Eqs. (jlHl), (ESI), and the dominant 
contribution to the surface integral in Eq. comes from Roafib- As expected, the phase 
factor depends on the Riemann curvature tensor, and is independent of coordinate choice, 
even though Eq. is dependent on the gauge- dependent field Na- 

We note also the underlying topological nature of Fig. Eb. If a path F is chosen that 
does not encircle the GW beam, AG = 0; if it does, AB 0. F cannot be shrunk to zero 
without cutting though the GW beam, and thus altering the topology of the F-GW beam 
system. Moreover, each time one goes around F, the phase difference AB changes by an 
integral multiple of the same factor; AB is thus proportional to the linking number of F 
around the GW beam. That AB is related to the linking number is not unexpected since Na 



is related to Ashtekar's 
experiments such as j^i 



oop variables. Nonetheless, it does point out the possibility of doing 
I that directly measures the linking number of F, and this indicates 



the underlying topological nature of this gravitational Aharonov-Bohm-type effect. 
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FIG. 7a FIG. 7b 

FIG. 7: Fig. 7a is a sketch of spacetime diagram showing the passage of the unipolar, Gaussian 
GW (shown in gray) and the worldhnes of two particles beams, along the two paths and 
(i. e. the two dominant Feynman paths), of the interferometer encircling it. Events pA and ps 
correspond to the splitting and recombining of the initial and final beam splitters shown in Fig. 7a. 
Fig. 7b shows a schematic of a possible experimental set-up (a Bonse-Hart interferometer) at the 
given time slice tj shown in Fig. 7a. The GW beam passes through the center of the Bonse-Hart 
interferometer, and the gray oval shows the cross section of the Gaussian beam waist which forms 
at the center of the plane of the interferometer. 

APPENDIX A: REVIEW OF LINEARIZED GRAVITY AND TETRAD FOR- 
MALISMS 

In this appendix, we present a brief review of some of the properties of general tetrad 
frames, differential forms, and linearize d g ravity that we shall need. The review is not 
exhaustive and the reader is referred to [2^ or for a more complete presentation. 

A tetrad is a local coordinate system formed by a set of four orthonormal vectors {e^^} 
such that Cj^^cb^ = tiab and ^cau = duu- These {e^'^} do not have to be tied to any 
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observer's worldline as we have done in the above, and the results presented here are vahd 
in general. 

Since tj^ is a unit vector, the covariant derivative V^e^,^ can only be a rotation or boost 
of e^^ . This "rotation matrix" is determined by the Ricci rotation coefficients: 

^CAB = ec^e^^V^CBu, (Al) 
with ucAB = —^CBA- It is straightforward to show that in a tetrad frame, 

RaBCD = e.j^^d^UJBCD — ^B^'^fJ.^ACD — 

{^AFC^B D ~ '^BFC^A D + '^AFB'^ CD ~ '^BFA^ CD ] ■ (A2) 

We emphasize that objects with capital Roman indices in the tetrad frame — such as 
^CAB — are scalars. They only take the partial derivative, and not the covariant deriva- 
tive. 

Given a tetrad frame, it is natural to work with differential forms, which we shall denote 
by symbols in boldface. As 1-forms we have = e^^dx^ and ujab = ^^^lABdx'^, and as 
2-forms we have Rab = R/^iuABdx^ A dx'^/2. As usual, d is the exterior derivative. Equation 
(IA1|) then becomes 

dsA = A ujab, (A3) 

where A is the wedge product; cj^^'s role as a "rotation" or boost matrix is now manifest. 
Equation ()A2jl is then 

R^^ = doj^^ + uj^^ Au)c^. (A4) 
Taking the exterior derivative of Eq. ()A3|) . we get the 1st Bianchi identity, 

= eij A Rj^^, (A5) 

and the exterior derivative of Eq. ()A5|) gives the 2nd Bianchi identity 

dR^^ = R/^ Au^c^ -uj/^Rc^. (A6) 

In the case of linearized gravity, g^i, = r]^y + h^y where h^^ is "small" , and we are only 
concerned with terms linear in /i^j,. Then g^'^ = rj^^ — h^^ where on the left hand side we 
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use rj^iy to raise and lower indices. In this limit 

R = -D/i^l + d^d.h^^". (A7) 

For GWs in the TT gauge, = 0, d'^h^^ = 0, and h^^ = 0, so that R = 0, Rq^ = 0, and 
the equation of motion for GWs is i?-^ = — = 0. 

To construct the tetrad frame for linearized gravity, we first note that in flat spacetime 
the tetrad frame is trivial: {5^}. The presence of small h^,^ rotates these vectors and 

= (^." - Ik) 5a. e^, = {^5; + \hi}j 5t. (A8) 

Note that this choice is not unique. As noted in Sec. IV 4, given any tetrad frame, we 
can always do a local Lorentz transformation that will still preserve the orthonormality of 
{e^'^}. The choice we have made for {e^^} for the GLF, and used in the construction in 
Sec. I, corresponds to an observer at rest in his proper frame. 

Using Eq. (jAlj) and the Levi-Civita connection in Eq. ()A7j) . we find that 

^CAB = ~ {dAliBc - dshAc) , (A9) 

while Rabcd = dA^scD — dsoJACD in the linearized gravity limit. Also in this limit, the 
equations for the Ricci tensor and scalar in Eq. ()A7j) have the same form in the tetrad frame, 
with the replacement of Greek indices by capital Roman indices. 



APPENDIX B: PHASE FACTOR SOLUTION 



Equation (jU^ is a quasi-linear partial differential equation [37j whose method of solution 
is well known. First, both t and X" are considered functions of a parameter r, so that 
9 = 6(t(r), X"(r)), and defines a constant surface in space. Consequently, 

dt dQ dX" ^ ^ 

and for Eq. ()B1|) to be a solution of Eq. (jHH), 

- 1 - AT- rR9l 

9^ - ^' ^ - ^ ' ^^^^ 
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the quasi-linear partial differential equation reduces to the solution of a set of ordinary 
differential equations. Although Eq. ()B2|) can be solved using standard methods once an 
initial condition is given, the underlying physics become much clearer if we consider instead 
the following function 

e~,{t{T),X'^iT)) = jl^^NMr),X'^iT))dX^, (B3) 

where the integral is from a fixed point {t{0), X°'{0)) to the point (t(T), X'^(T)) along F. 
The prefactor m/h is included so that Qp is unitless. We choose F such that its tangent 
vector is given by Eq. ()B2|) . and we parameterize it by r G [0,T]. F starts at the point 
(t(0), X"(0)), and is used also as the initial condition for Eq. ()B2j) . Clearly, 



def{t,X'') _m dX 



A 



m 



dr h dr h 

in the linearized gravity limit, while the left-hand-side vanishes identically from Eq. ()B1|1 . 
Thus, Qf{t,X"-) is a solution of Eq. (jHSj)- Because it is path dependent, this solution is not 
unique, however. 

The spatial component of the tangent vector to F lies along A''^, and for plane waves, Na 
is perpendicular to the direction of propagation. Thus, for the Gaussian beam in Fig. [7|3, F 
will wrap around the beam. We can, of course, go around this beam either in a clockwise or 
counterclockwise direction, and we denote a clockwise path by F^ and a counterclockwise 
path by F^ . To F~ there is then the corresponding the function 0j,_, and to F^ there is 
the corresponding function Qf+. 

Consider now the functions 

e^i(t(T),X«(T)) = "^j^NAdX^, 

eMt{T),X^{F)) = ^ / NAdX^, (B5) 

h jp2 

and we restrict ourselves to those F^ that lie on the surface of the beam, meaning that F^ 
do not come closer than lew to the center of the beam of GWs. Because Na is exponentially 
small outside of the Gaussian beam, G^i ~ Qp+ and 9/^2 ^ Qp_. Thus, Eqs. flB5|) are 
solutions of Eq. (jU^ . 
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